clc,clf,clear;
f=chebfun('1./(1+x.^2)',[-5,5]);
plot(f);
p=chebpoly(f);